Ferromagnetism in Carbon doped Zinc Oxide Systems 
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We report spin polarized density functional calculations of ferromagnetic properties of a series of 
ZnO clusters and solid containing one or two substitutional carbon impurities. We analyze the eigen 
value spectra, spin densities and molecular orbitals, and induced magnetic moments for ZnC, Zn2C, 
Zn20C, carbon substituted clusters Zn„O n (n=3-10, 12) and ZnO solid. The results show that the 
doping induces magnetic moment of the ~2 /is in all the cases. All the systems with two carbon 
impurities show ferromagnetic interaction except when the carbon atoms share the same Zn atom 
as the nearest neighbor. This ferromagnetic interaction is predominantly mediated via n bonds in 
ring structures and through it and a bonds in three dimensional structure. The calculations also 
show that the interaction is significantly enhanced in solid, bringing out the role of dimensionality 
of Zn-O network connecting two carbon atoms. 

PACS numbers: 36.40.-c, 71. 15. Mb, 75.50.Pp 



I. INTRODUCTION 

Dilute magnetic semiconductors (DMS) are a focus 
of much attention due to their potential application as 
spintronics material. These materials, apart from hav- 
ing a desirable band gap can exhibit room temperature 
ferromagnetism (RTF) upon doping by suitable impuri- 
ties. Several oxides like TiC>2, SnC>2, In203, Hf02 and 
more popular ZnO have been investigated and have been 
found to show ferromagnetism when doped with transi- 
tion metal impurities Amongst these, ZnO is pre- 
ferred because of many attractive features. It has a wide 
band gap of 3.37 eV, suitable for applications with short 
wavelength light, a large exciton binding energy, a low 
lasing threshold, and is transparent to visible light. 

There are a number of theoretical and experimental 
studies reported on ZnO yielding a variety of data on 
magnetic behavior. The data has raised many issues and 
at times, controversies. The density functional calcula- 
tions show that doping by Cu produces half metallic fer- 
romagnetpi, the prediction has been verified experimen- 
tally. 7 Co doping has generated much controversy about 
the origin of magnetism. Cobalt doped ZnO is reported 
to have high Curie temperature (TC) of up to 700 K. 
A robust ferromagnetic state is predicted at the Oxygen 
surface even in the absence of magnetic atoms4>- 9 If Mn 
or Co is doped in films of ZnO, ferromagnetic behavior 
has been observed in both insulating and metallic films, 
but not when the carrier density is intermediated 

Recently, non magnetic impurities have been shown to 
induce ferromagnetic behavior. Vacancy-induced mag- 
netism in ZnO thin films and nano wires have been re- 
ported on the basis of density functional investigations. 11 
In fact, defects of several types like vacancy at Zn site, 
oxygen surface, interstitial edges in nano ribbons** etc. 
have been reported to be inducing magnetic interactions. 
Recently, carbon doped ZnO has been shown to exhibit 



ferromagnetic behavior with Curie temperature around 
room temperature^ The experiment has been performed 
with two levels of doping, 1 % and 2.5 %. The spin 
density functional theory (SDFT) results show that the 
magnetism is due to ZnC system. The magnetic moment 
per carbon atom decreases slightly for the higher con- 
centration from 2-3 /is to 1.5-2.5 /is. There are several 
unclear issues: the origin of magnetic moment on carbon 
site, the role of Zn and oxygen, the nature of the long 
range interaction between the two carbon sites, the role 
of oxygen-Zn network connecting the two sites. In order 
to gain some insight in to some of these issues, we have 
carried out a systematic density functional investigations 
on a series of systems, beginning with ZnC molecule. We 
have investigated the evolution of magnetic moments and 
magnetic interaction in ZnC, Zn2C, Zn20C, Zn2C2 and 
then in a series of clusters Zn„O n (n=3— 12) with a single 
and double carbon substitution. Finally we examine ZnO 
solid with a large unit cell of 72 atoms. We bring out the 
role of exchange splitting in Carbon atom, hybridization 
of Carbon p orbitals with Zn, the role of oxygen-Zn net- 
work, especially the role of dimensionality in stabilizing 
magnetic interactions. 



II. COMPUTATIONAL METHOD 

The calculations have been performed using SDFT 
with ultra soft pseudopotential and plane wave basisi 14 ! 15 
We have used two different codes namely CASTEP 15 
and VASPd In all the cases, we have used generalized 
gradient approximation (GGA) given by Perdew-Burke- 
Ernzerhof (PBE)d For clusters we have used supercells 
having lengths ~15-30 A. The energy cut-off for the 
plane wave basis was chosen to 300 eV. The structures 
were considered to be converged when the force on each 
ion was less than ~0.05 eV/A with a maximum dis- 
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FIG. 1: Optimized ground state geometries of C-doped (ZnO) n (n=l-5) clusters. The binding energy per atom (E&) in eV is 
also shown. The oxygen atoms are indicated by red color (online) whereas substituted carbon is marked by C. The second and 
third column shows structures with single and double carbon substitution respectively. 



placement of ^0.002 A and a convergence in the total 
energy of about 2 x 10~ 6 eV/atom. In all the cases, the 
geometries of all the lowest energy structure are in good 
agreement with results of Matxain et. al. 18 To study 
the magnetic properties of C-doped (ZnO)„ clusters, car- 
bon atoms were doped substitutionally at various O-sites. 
The C-doped clusters were also optimized with the same 
set of parameters as used for (ZnO)„ clusters. 



We have verified the accuracy of our model on ZnO 
molecule and ZnO solid. We find a bond length of 
1.736 A and a binding energy 1.79 eV for ZnO molecule. 
The previous studies^ also report a bond length between 
1.71-1.75 A. We also calculate the binding energy of the 
ZnO solid in the wurtzite structure. The calculated bind- 
ing energy of 7.23 eV per ZnO unit agrees well with the 
experimental binding energy of 7.52 eV per ZnO as re- 
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FIG. 2: Optimized ground state geometries of C-doped (ZnO) n (n=6-8, 12) clusters. The binding energy per atom (Ej,) in eV 
is also shown. The oxygen atoms are indicated by red color (online) whereas substituted carbon is marked by C. The second 
and third column shows structures with single and double carbon substitution respectively. 



ported by Jaffe et.al^S 

In order to understand the basic mechanism of ferro- 
magnetic origin and evolutionary trends, we have per- 
formed the calculations on series of clusters namely; 
ZnC, Z112C, (ZnO)„ and the substitutionary doped 
Zn„0„_, n C m clusters with n = 2 — 10, 12 and m = 1 — 2 
as well as on ZnO solid. In each of the cases, we substi- 
tute either one or two carbon atoms at O-sites. 



In next section, we present and discuss our results. 



III. RESULTS AND DISCUSSIONS 

It is most fruitful to examine the electronic structure 
and systematics of magnetic properties of some impor- 
tant motifs in zinc oxide solid. When the carbon atom is 
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FIG. 3: The Eigenvalue spectra for some clusters. In each 
of the cases the left column is for up and right column is for 
down electrons. The unoccupied levels are shown by dotted 
lines. The numbers near the line shows degeneracy. 



substitutionally doped at O-site in the bulk zinc oxide, 
it has zinc atom as the first nearest neighbor and oxygen 
atom as the second nearest neighbor. Therefore, we be- 
gin by examining the magnetism and the electronic struc- 
ture of ZnC, Zn2C, Z^OC, and Zn2C2 clusters. We have 
also calculated the electronic structure and the magnetic 
properties of ring structures of Zn„O n cluster (n=3-8) 
and ZngOg-Zni20i2 which form three dimensional struc- 
tures. The results are also presented for bulk zinc oxide 
by using large unit cell with 72 atoms. The optimized 
ground state geometries of these ZnO systems are shown 
in figure 1 and 2. It can be seen that the substitutional 
doping by carbon atoms has a small effect on the geome- 
tries of the hosts with exception of Zn2 0. 

In table 1, we summerize the results of all C-doped 
clusters. The table shows the magnetic moments per site, 
the exchange energy Ej calculated as Efm-^afm (wher- 
ever appropriate) where, ~Efm an d ^afm are the total 
energies of ferromagnetic and antiferromagnetic state re- 
spectively, the nature of magnetic ground state and na- 
ture of ground state geometries of C-doped clusters. A 
number of features can be discerned from table 1. 1. A 
majority of doubly doped (ZnO)„ clusters turn out to be 
ferromagnetic with magnetic moment of ^2 fiB per car- 
bon atom. 2. Whenever the same zinc atom is the nearest 
neigbhor to both the carbon atoms, the ground state is 
antiferromagnetic. This is the case for Zn2C2, Zn30C2 
and in other clusters when the carbon atoms are substi- 
tuted on both sides of zinc atom. 3. It turns out that 
^80-90 % magnetic moment is localized on carbon sites 
with small but significant induced magnetic moment on 
oxygen and zinc sites (see figure 5 and the discussion). 
This is true for all singly and doubly substituted clusters 
irrespective of the placement of carbon atoms. 4. A re- 



markable observation is that the exchange energy (Ej) is 
sensitive to the dimensionality of clusters and is highest 
for three dimensional clusters. Thus, the number of paths 
connecting two carbon atoms via zinc-oxygen network is 
a significant factor favouring the ferromagnetic coupling. 
In a sense this also brings out the role of oxygen in es- 
tablishing the ferromagnetic interaction. We elaborate 
on this point in the discussion below. 

In order to elucidate the origin of magnetism, we ex- 
amine the eigenvalue spectra (figure 3) and the nature 
of molecular orbitals for ZnC, Zn2C, Zn20C and Zn2C2 
along with atomic Zn, O, and C which are shown in fig- 
ure 4. We use the nomenclature alpha homo and beta 
homo to denote the highest occupied molecular orbitals 
for up spin and down spin respectively. In figure 3, we 
show the eigenvalue spectra for up spin (left part) and 
down spin (right part) electrons. The relevant molecular 
orbitals are shown in figure 4 as isovalued surfaces. We 
also show the induced magnetic moments on all sites in 
figure 5 for some of the relevant systems. We will dis- 
cuss the spectra, molecular orbitals and site distribution 
of the magnetic moments together. Let us note that the 
exchange spliting for carbon atom between spin up and 
spin down electron is ^2.49 eV. The magnetic moment is 
due to two up spin occupied orbitals, consistent with the 
Hund's rule. This brings the eigen values of atomic car- 
bon p and zinc 4s quite close to each other for up states, 
while the down spin eigen value of carbon p is placed 
significantly higher. The formation of the spectra across 
the series reveals that the action is in the s-p complex 
formed by carbon p and zinc 4s orbitals, and the spliting 
of the up and down "bands" persists for all the systems, 
leading to the formation of the magnetic moment around 
carbon. 

ZnC carries a magnetic moment of ~2 /ib on carbon 
site. Because of the linear nature of ZnC, only one of the 
p orbitals is capable of hybridizing with Zn s. As a conse- 
quence triply degenerate p states (for the both the spins) 
split into two states as a singly degenerate and doubly de- 
generate once, retaining the exchange splitting between 
up and down. The corresponding molecular orbitals are 
shown in figure 4. A lowest state in the complex for 
both the spins are of similiar nature and show significant 
hybridization and delocalization. The magentic moment 
is due to the doubly degenerate pure p states localized 
on the carbon atom. The Mulliken population analysis 
shows effective charge transfer from Zn to C of order of 
~1.7 electrons. Thus ignoring the splitting between up- 
down eigenvalues, a simple picture emerges. The charge 
transfer from zinc makes one of p orbitals on carbon dou- 
bly occupied. The delocalization caused by hydridization 
with 4s reduces Coulomb repulsion. The magentic mo- 
ment is due to the doubly degenerate pure p states local- 
ized on the carbon atom. 

Next we turn to Zn2C. Since in ZnC, two spin down 
orbitals on carbon are still unoccupied, one would naively 
expect Zn2C to be nonmagnetic. Interestingly, Zn2C 
also shows a magnetic moment of ~2/is mainly local- 
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TABLE I: The magnetic moment (fiia) in /is, Exchange energy Ej (Energy difference between ferromagnetic and antiferro- 
magnetic state) in meV, the nature of magnetic ground state for all the systems investigated. 



Host Systems 


Magnetic Moment jiM 


Ej (meV) 


Nature of 
prmiTid state 


Geometry 


Singly 
doped 


Doubly 
doped 




2 


- 








ZnO 


2 


- 






Linear 


Zn2<J2 


2 







At* 1V1 


Rhomb ous 


ZmOs 


2 







AFM 


Cap-triangle 


ZJYI4U4 


2 


4 


-26.4 


1 M 


Square 


ZnrjOr, 


2 


4 


-3.3 


1 M 


Ring 


ZneOe 


2 


4 


-2.5 


1 M 


Ring 


Zn 7 7 


2 


4 


-1.1 


FM 


Ring 


Zn 7 7 a 


2 






AFM 


Ring 


Zn 8 O s 


2 


4 


-0.9 


FM 


Two connected rings 


Zn 9 9 


2 


4 


-46.3 


FM 


3D 


ZnioOio 


2 


4 


-26.3 


FM 


3D 


Zni 2 Oi2 


2 


4 


-85.3 


FM 


3D 


ZnO soild b 


2 


4 


-89.3 


FM 


Wurtzite 



"Zinc is nearest neighbour to both carbon atoms. 
'Calculation was performed on ZnO solid within large unit cell 
containing 72 atoms. 



ized around the carbon atom. The hybridization of one 
of the p orbitals with two zinc 4s splits the spectrum of 
up and down electrons forming bonding and antibonding 
orbitals, the splitting being much less in the case of down 
electrons. As a consequence four up and two down or- 
bitals are occupied. Once again by examining the molec- 
ular orbitals (figure 4) , we arrived at a simple description 
namely; there are two doubly occupied orbitals which are 
delocalized (2p-4s hybridized). The magnetic moment 
arises out of two carbon based p's orbitals (homo-1 and 
homo-2 for up electrons). As a consequence of the hy- 
bridization a small but significant magnetic moment is 
induced on both zinc atoms. 

Next we examine a typical motif Zn 2 OC which brings 
in the role of oxygen. With introduction of oxygen, there 
is a competition between oxygen and carbon for charge 
transfer from zinc. The occupied molecular orbitals show 
interesting pattern. Evidently, relevant orbitals leading 
to magnetic interactions are homo and homo-1. The 
homo shows tt bond formation between carbon and oxy- 
gen (homo perpendicular to plane) while homo-1 shows 
the a bond. There are three down orbitals nearly com- 
pensating the charge in the three up orbitals. It can be 
seen that a significant magnetic moment (^0.42 /ig) is 
induced on the oxygen sites, clearly because of homo and 
homo-1. The existence of 7r and a bonds establishes the 
coupling between carbon and oxygen and is a crucial ele- 
ment in the formation of long range magnetic interactions 
even in the extended systems. 

The ground state of symmetric cluster Zn 2 C2 is anti- 
ferromagnetic. There is a charge transfer of one electron 



to carbon atoms from each of the zinc atoms. As a conse- 
quence the zinc-zinc distance increases by ^0.2 A. As can 
be seen from figure 4, it is homo-2 which is responsible 
for antifcrromagnctic coupling. This is analog of super- 
exchange mechanism. It may be noted that a similiar 
observation can be made for Zn 3 OC2 as well as for all 
the ring clusters when carbon atoms are placed on both 
sides of zinc atom. We recall that the experiment indi- 
cate a saturation of magnetic moment per carbon atom 
around 5 % substitution. 

Now, we turn our attentation to ring structures Zn 5 C>5 
- Zn 7 7 . We discuss the typical case of Zn 7 05C2 show- 
ing the ferromagnetic behaviour. The relevant molecular 
orbitals are shown in figure 6. In this case two carbon 
atoms are separated by chain of zinc-oxygen-zinc (sec 
figure 5). It can be seen from alpha homo and alpha 
homo-1, that there is a ir bond formation between three 
p orbitals centered on two carbon atoms and one oxygen 
atom. A similiar bonding is also seen in alpha homo-2 
and alpha homo-3 as well as in beta orbitals which are 
shown in figure 6. These are the orbitals which are re- 
sponsible for the ferromagnetic interaction between the 
carbon atoms mediated via oxygen. 

Zni20ioC2 is a three dimensional structure. As 
pointed out earlier when the structure changes to three 
dimensional one, there is a significant change in the con- 
nectivity of two carbon atoms via zinc and oxygen. The 
ground state structure for this system is shown in figure 
2. The induced magnetic moments are shown in figure 5. 
It can be seen that the induced magnetic moments on two 
nearest neighbor oxygen sites are substancial ^0.180 fis- 
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FIG. 4: Isovalued surfaces of molecular orbitals for up(alpha) and down(beta) electrons. The position of substituted carbon is 
marked by C. 



Other oxygen sites connecting them also have notable 
magnetic moments. The relavant molecular orbitals (al- 
pha homo - alpha homo-2 ) and the spin density are 
shown in figure 7. Clearly the picture described in the 
earlier cases continues in the present case. The induced 
spin density shows the magnetic moment on carbon and 
oxygen sites. These homos show interaction between two 



carbon atoms mediated via oxygen. Thus, in these case 
once again the magnetic interaction between the carbon 
atoms is mediated via p orbitals (7r-bond). The most 
significant point is that the interaction is strenghtened 
by availability of multiple sites connecting two carbon 
atoms. In this case, the C-C distance is ~ 5.410 A . 

We expect a similiar picture to be valid in ZnO solid. 
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FIG. 5: Spin moment at site is depicted in C-doped clusters. The position of substituted carbon is marked as C. 



We have considered a large unit cell with 72 atoms con- 
taining two carbon atoms separated by 7.6 A. Our results 
for the density of states (DOS) and magnetic moments 
are consistent with the report of of Pan et.alM- We find 
that the spin-up bands are fully occupied while the spin- 
down bands are partially filled resulting in the magnetic 
moment of ^2/is per carbon atom. We show the relavant 
molecular orbitals near the top of the valence band for 
up electrons. It can be noted from figure 8 that the two 
carbon atoms are separated by a zinc oxide plane. We 
also show the spin density in figure 9. The molecular or- 
bitals depict two different ways in which the interaction 
is mediated. In the case of alpha- 1, the carbon (in the top 
plane) induces magnetic moments on oxygen sites within 
the plane which then couples to the oxygen in the bottom 
plane as seen in figure 8. In the case of orbitals marked 



as alpha-2, the carbon induces the magnetic moments 
in the plane, below and above which couples the carbon 
site in the bottom plane. These figure clearly bring out 
the role of three dimensional oxygen network in provid- 
ing the multiple paths between two carbon atoms. The 
resulting spin density is shown in figure 9. As expected 
the spin density is positive and is mainly around carbon 
sites. The induced magnetic moments on oxygen sites 
are clearly seen. 



IV. SUMMARY AND CONCLUDING 
REMARKS 

We have carried out spin polarized electronic structure 
calculations on series of (ZnO) n and the substitutionally 
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FIG. 6: Isovalued surfaces of molecular orbitals for up and down elecrons in C-doped ZnyOy clusters. The oxygen atoms are 
indicated by red color (online) whereas the position of substituted carbon is marked by C. 



doped Zn n O n - m C m clusters with n = 1 — 10, 12 and 
m = 1 — 2 as well as for zinc oxide solid with a view to 
elucidate the origin of magnetism in carbon doped zinc 
oxide. Our calculations reveal the crucial role played by 
oxygen in mediating ferromagnetic interaction. A de- 
tailed study of molecular orbitals show that the inter- 



action is mediated via it as well a bonds (especially in 
three dimensional solid) formed out of p electrons of car- 
bon and oxygen. Our calculations also bring out the role 
played by hybridization of zinc 4s with carbon 2p and 
oxygen 2p. We expect the antiferromagnetic coupling 
if carbon atoms get substituted on sites adjacent to zinc 





Alpha homo-2 Spin density at site 



FIG. 7: Isovalued surfaces of molecular orbitals for up and down electrons in C-doped Zni20i2 clusters. The isosurface spin 
density at site is also shown. The position of substituted carbon is marked by C. 



leading to saturation of magnetic moment with increasing 
carbon content. The results clearly show a significant en- 
hancement in the ferromagnetic interaction between the 
two carbon atoms due to higher dimensionality, bringing 
out the role of multiple connectivity (between two carbon 
atoms) of the Zn-0 network. 
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